Double resonance in the infinite-range quantum Ising model 
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We study quantum resonance behavior of the infinite-range kinetic Ising model at zero temper- 
ature. Numerical integration of the time-dependent Schrodinger equation in the presence of an 
external magnetic field in the z direction is performed at various transverse field strengths g. It is 
revealed that two resonance peaks occur when the energy gap matches the external driving frequency 
at two distinct values of g, one below and the other above the quantum phase transition. From the 
similar observations already made in classical systems with phase transitions, we propose that the 
double resonance peaks should be a generic feature of continuous transitions, for both quantum and 
classical many-body systems. 
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I. INTRODUCTION 

A noise is often considered as a nuisance for a system 
to display any ordered behavior, and thus the weaker 
the better for the performance of the system. However, 
for the last decades, a lot of researchers have revealed 
that this is not always the case and that there exist a 
class of systems in which the intermediate strength of 
noise can help the system to show the best coherence 
with an external periodic driving. This surprising phe- 
nomenon was termed as the stochastic resonance (SR) 
due to its stochastic nature The phenomenon of SR 
has been found in the fields of physics and earth sci- 
ence, as well as in biology: the periodically recurrent 
ice ages, the bistable ring laser, superconducting quan- 
tum interference device, human vision and th3 auditory 
system, and the feeding mechanism of paddle fish, to 
list a few [l|, The occurrence of the SR is properly 
explained by the time-scale matching condition: the co- 
herence between the system's response and the external 
driving becomes strongest when the stochastic time scale 
inherent in the system matches the time scale provided 
by the external driving. In a simple classical system of 
a few degrees of freedom making contact with a thermal 
reservoir, the intrinsic time scale is given by the mono- 
tonically decreasing function of the exponential thermal 
activation form. Accordingly, the above mentioned time- 
scale matching condition can only be satisfied at a single 
temperature [lj. The time-scale matching condition was 
later extended to the classical statistical mechanical sys- 
tems with continuous phase transitions such as the glob- 
ally coupled, i.e., infinite-range kinetic Ising model Q. 
It has been shown that the nonmonotonic behavior of 
the intrinsic time scale around the critical temperature 
makes the time-scale matching condition satisfied at two 
distinct temperatures, one below and the other above the 
critical temperature, resulting in the double resonance 
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FIG. 1. Infinite-range (a)-(b) classical and (c)-(d) quantum 
Ising models. Double SR peaks are observed in (a) the oc- 
cupancy ratio (OR, the fraction of spins in the direction of 
the external periodic driving), and in (c) the power spectrum 
at the driving frequency Q, Pn(g) = \J dte' Ljt m(t; g)\u=ni 
where m(t; g) is the magnetization. Fluctuation strength is 
controlled by (a) the temperature T and (c) the transverse 
field strength g, respectively. Arrows denote the positions of 
maxima in (a) and (c), which are in good agreements with the 
time scale matching conditions shown in (b) and (d). In (b) r 
is the intrinsic relaxation time scale, and in (d) A is the energy 
gap (and thus the inverse time scale). The dotted horizontal 
lines in (b) and (d) denote the external frequency scales of 
time-periodic driving. A in (d) is obtained from the direct 
diagonalization of the quantum Hamiltonian for the system 
size N — 1000 and Pn (g) in (c) through the semiclassical ap- 
proximation on the Heisenberg equation of motion (see text). 
See p| for (a) and (b). 



peaks. The double SR peaks have also been observed in 
the classical Heisenberg spin system in a planar thin film 
geometry and infinite-range g-state clock model [Bj]. 

The SR phenomenon in quantum systems, named as 
quantum stochastic resonance (QSR), has been studied 
with focus on the interplay between quantum and clas- 
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sical fluctuation at finite temperatures 0, Q. The QSR 
at zero temperature has also been studied for the one- 
dimensional quantum spin system with a spatially modu- 
lated external field, and the length-scale matching similar 
to the time-scale matching in conventional SR has been 
discussed pj. In the present work, we study the QSR 
at zero temperature in the Ising spin system with the 
quantum phase transition [8j. We summarize our main 
findings in Fig. [T] which displays the double SR peaks 
and the time-scale matching conditions in infinite-range 
classical [3| and quantum (this work) Ising systems in the 
presence of a weak external driving with the frequency 
O. We conclude that the time-scale matching condition 
allows us to understand the classical and the quantum 
double SR peaks on the same ground. 

In this paper, we numerically study the resonance be- 
havior of the infinite-range quantum Ising model. Inte- 
grations of the time-dependent Schrodinger equation and 
the semiclassical equation of motion unanimously yield 
the existence of the double SR peaks, which are clearly 
explained from the matching condition between the en- 
ergy gap, which is intrinsic, and the frequency of the 
external time periodic driving. 
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FIG. 2. (Color online) The energy gap A between the ground 
and the first-excited states versus the transverse field strength 
g for the globally-coupled quantum Ising model without exter- 
nal driving in z-direction. A vanishes as the quantum critical 
point g c = 1 is approached. Inset: Finite-size scaling of the 
form A = {g - g c ) zv f[(g - g c )N^ d ^] with z = 1, v = 1/2, 
and d c = 3. 



II. RESULTS 

Let us begin with the globally coupled N 
spins described by the Hamiltonian H = 
-l/(2NS)'£# k S*Si - ;/V ; .s;. where Sf is the 
spin-f/2 operator in the a direction (a = x,y,z) at 
the jth site (S = 1/2 and h = 1 henceforth), and the 
transverse field g in the x direction induces quantum 
fluctuation due to [Sj,S^] = iSjkS^ ^ 0. By using the 
total spin operator J a = J^j S° with J = N/2, the 
Hamiltonian can be cast into the form [9fl 
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which allows us to handle much bigger N since the num- 
ber of base kets becomes only N + 1 (we use J z eigenkets 
as base kets) . The globally coupled quantum Ising model 
Eq. ([TJ is very well-known to exhibit the quantum phase 
transition of the mean-field nature and its finite-size scal- 
ing has also been extensively studied @ (see [l(| for the 
finite-size scaling of the quantum phase transition in the 
one-dimensional Ising chain system). 

We numerically obtain the energy gap A between the 
ground and the first-excited states of the Hamiltonian 
Eq. (UJ, which exhibits the quantum phase transition at 
g = g c = 1 of the mean-field universality class as dis- 
played in Fig. [2] for the system sizes N = 200, 600, and 
fOOO. The inset of Fig. [2] shows the finite-size scaling of 
A with the well-known exponents: dynamic critical ex- 
ponent z — I, the correlation length exponent v = 1/2, 
and the upper critical dimension d c = 3 Q. The vanish- 
ing energy gap (and thus the divergence of the intrinsic 
time scale) at the quantum critical point is particularly 



important in the present study: The nonmonotonicity of 
A as a function of the fluctuation strength g provides 
the origin of the double quantum resonance peaks (see 
Fig. P. 

In parallel to studies of the classical SR behaviors 
we next apply the weak time-periodic external magnetic 
field hit) — ho cos fit along the z-direction with ho = 
10~ 3 and £1 = 0.8, to get the time-dependent Schrodinger 
equation 
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dt 



N 



J 2 Z - g,J x - h(t)J z 



!*(*)>, (2) 



where the quantum ket l'F(i)) = Elf=-j^w(t)|-M] 
with J Z \M) = M\M) and the complex coefficient 
Ajvf(t). The time evolution of the system is numerically 
traced through the use of the fifth-order Runge-Kutta 
method combined with the Richardson extrapolation and 
Bulirsch-Stoer method [ll|. We check that the use of 
the sufficiently small time step 5t = 10~ 4 keeps the nor- 
malization condition |^4m(0| 2 = 1 unchanged within 
numerical accuracy. We first get the ground state in the 
presence of the extremely small external field in the pos- 
itive z-direction to break the up-down spin symmetry, 
and use it as the initial condition for Eq. ©. 

As the most important quantity to detect SR behav- 
ior, the average magnetization in the z direction m(t) = 
(1/J){^j (t)\ J z \^>(t)) is measured as a function of time. 
We do not observe significant difference for other sys- 
tem sizes, and we display our results m(t) for N = 600 
in Fig. [3] at (a) below and (b) above the quantum crit- 
ical point g c = 1. When g > g c , m(t) oscillates around 
m = 0, and we shift vertically each m(t) in Fig. [H^b) 
for better comparison. It is obvious from Fig. [3] that 
the resonance behavior of m(t) is seen in the form of 
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FIG. 3. (Color online) The order parameter 

m(t)=(l/iV)<*(t)|J*|*(t)) in time t at (a) g < g c = 1 and 
(b) g > g c for iV=600 in the presence of an external oscillat- 
ing magnetic field ho cos Qt with S7=0.8 in the z direction. At 
two values of g, one at gi ~ 0.59 and the other at <?2 — 1.44, 
m(t) exhibits the larger oscillation amplitude reflecting the 
stronger coherence with the external driving. For better dis- 
tinction, m(t) is vertically shifted by 0.9, 0.7, 0.4, and at 
g = 1.10, 1.40, 1.44, and 1.50 in (b), respectively, (c) The os- 
cillation amplitude Am versus g, exhibiting double SR peaks. 



the larger oscillation amplitude at two distinct strengths 
of quantum fluctuation, i.e., one below g c and the other 
above g c . In Fig. (3Jc), we display the oscillation ampli- 
tude Am = max t m(t) — min t m(i), which clearly shows 
double resonance peaks. We denote the first and the sec- 
ond resonance points as g\ « 0.59 and g 2 ~ 1-44, where 
the oscillation amplitudes become maxima. As another 
indicator of the SR behavior, we carry out the Fourier 
transformation of m(t) to obtain m{uj) at frequency uj. 
Figure |4] displays the magnitude of the spectral compo- 
nents 1 771 (w) I versus u) at various values of g. In general, 
there exist two peaks in |m(w)|, one at u)% = Q = 0.8 
[indicated by the dotted vertical lines in Fig.0], and the 
other at the position uj 2 that depends on g. We observe 
that the latter peak at lo 2 simply originates from the en- 
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FIG. 4. (Color online) The magnitude of the power spectral 
density |m(u)| versus frequency u) at various values of the 
transverse field strength g. The dashed vertical lines denote 
the frequency of the external field (Q. = 0.8). At g = 0.59 and 
1.44, the two peaks in the spectral density merge to become 
one. The peak in (b) at ui=l.Q is a harmonics, and does not 
reflect the energy gap like peaks in other panels. 



ergy gap (see Fig. [5]), i.e., uj 2 = A (note that H = 1 in 
this work). As g is increased toward g c from below, A 
decreases (see Fig. [2]), which in turn yields decreasing u> 2 
as shown in Fig. H{a)-(c). As g passes through g c from 
below, L02 bounces back and moves to a larger value, re- 
flecting the increase of A for g > g c in Fig. [5] Another 
important observation one can make from Fig. [4] is that 
when the two peaks at u>i and 0J2 merge into a single one 
at cj = uj\ = Q, the power spectrum (P w = |m(w)| 2 ) at 
fl suddenly increases much. From this, it is clear that 
the merging of the two peaks must occur at two distinct 
values of g, which are in good agreement with g\ ss 0.59 
and g 2 f» 1.44 in Fig. [3] 

We next adopt the Heisenberg picture in which the 
spin operator satisfies the equation of motion J a — 
— i[H, J a \. By using the commutation relation [J a , Jp] = 
ie a p 1 J 1 , we get J x = (1/2 J)(J z J y + J y J z ) + h(t)J y , J y = 
-(1/2J)(J Z J X + J X J Z ) + gj z - h(t)J x , and j z = -gJ y . 
We then make the semiclassical approximation and treat 
the spin operator J a as the ev-th component of the classi- 
cal spin vector J = J(sin 6* cos c6, sint9sino/>, cost9), which 
results in @, QJ] 



g sin <p. 

g cot 9 cos <t> — cos 8 — h(t). 



(3) 



We take initial values of 9 and <fi from the ground state 
of the system calculated by semiclassical approximation 
of Eq. ([T]) , and perform numerical integrations of Eq. Q 
in time at given values of g. In this semiclassical ap- 
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(a) g < g c 




FIG. 5. (Color online) |m(cj)| versus g and w in the two re- 
gions (a) g < g c = 1 and (b) g > g c — 1 obtained via semiclas- 
sical approximation applied for the Heisenberg equations of 
motion. The right panels of (a) and (b) show the color scheme 
for |m(w)|. Along uj — 0.8 axis, as g is increased, |m(w)| ex- 
hibits the biggest value both at g ~ 0.59 and g ~ 1.44. 



two distinct values of g: g\ ~ 0.59 and gi ~ 1.44 as dis- 
played in Fig. [5j which are in perfect agreement with the 
findings in Fig. 2) 



III. SUMMARY 

In summary, the infinite-range quantum transverse- 
field Ising model at zero temperature has been numer- 
ically investigated in the presence of the weak longitudi- 
nal time-periodic magnetic field at the frequency fi. The 
resonance behavior at two distinct values of the trans- 
verse field g has been clearly observed via (i) the large 
amplitude oscillation of the magnetization in time and 
(ii) the large peak at Q in spectral analysis. The origin 
of the double SR peaks in the system has been identified 
from the vanishing of the energy gap around the quan- 
tum critical point. When the energy gap matches the fre- 
quency of the external field, the strong resonance peaks 
occur at two different values of <?, exhibiting the double 
resonance behavior. We have also confirmed the double 
resonance in the thermodynamic limit through the use 
of the semiclassical approximation made for the Heisen- 
berg equation of motion. We propose that the time-scale 
matching condition should play an important role in un- 
derstanding the double SR behavior in a broad range of 
systems with continuous phase transitions, both classical 
and quantum. 
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